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ABSTRACT 

We investigate the hadron collider phenomenology of the Minimal Supersymmetric Standard 
Model (MSSM) with explicit CP violation for Higgs bosons that can be observed in Standard 
Model search channels: W/hh) at the Tevatron, and gg —> Hi{-^ 77 ), ttHi{-^ bb) and 
WW Mt~) at the LHC. Our numerical analysis is based on a benchmark scenario 

proposed earlier called CPX, which has been designed to showcase the effects of CP violation 
in the MSSM, and on several variant benchmarks. In most of the CPX parameter space, these 
hadron colliders will find one of the neutral MSSM Higgs bosons. However, there are small regions 
of parameter space in which none of the neutral Higgs bosons can be detected in the standard 
channels at the Tevatron and the LHC. This occurs because the neutral Higgs boson with the 
largest coupling to W and Z bosons decays predominantly into either two lighter Higgs bosons 
or a Higgs boson and a gauge boson, whilst the lighter Higgs boson has only small couplings to 
the W and Z bosons and the top quark. For other choices of CP-violating parameters, all three 
neutral Higgs bosons can have significant couplings to W and Z bosons, producing overlapping 
signatures: these may or may not be distinguishable from backgrounds. The existence of these 
regions of parameters provides a strong motivation for a detailed experimental simulation of these 
channels. 
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1 Introduction 


Detailed studies of the mechanism of electroweak symmetry breaking top the agendas of 
operating particle accelerators, as well as those under construction or being proposed for 
the near future. One paradigm for many of these studies is provided by the minimal super- 
symmetric extension of the Standard Model (MSSM) There is increasing interest in the 
possibility that the MSSM includes explicit CP violation, which may provide opportunities 
to probe CP violating parameters through the Higgs sector 


m, 11,12,13 


Assuming universality of the soft supersymmetry (SUSY)-breaking gaugino masses 
mi/ 2 , scalar masses mo and trilinear parameters Aq at the Grand Unihcation Theory (GUT) 
input scale, the MSSM then provides only two new sources of CP violation in addition to 
the Kobayashi-Maskawa phase of the Standard Model (SM). In a suitable convention, these 
may be represented as the phases of complex parameters mi /2 and Aq. The values of these 
phases are constrained by the upper limits on the electric dipole moments (EDMs) of the 
electron, neutron |T^,|^ and mercury atom [Q and other measurements. The ongoing 
probes of CP violation at B factories can also test the consistency of the above minimal 


CP-violating supersymmetric scenario [17,18 


In the MSSM with explicit CP violation, the effect of CP violation on the Higgs sector 
enters beyond the Born approximation @,|^. The effects induced by arg(Ao) and 

arg(mi/ 2 ) on the MSSM Higgs-boson masses have been studied in some detail [^,^. Several 
interesting phenomenological implications for the Higgs bosons emerge in this minimal CP- 
violating framework of the MSSM. The three neutral Higgs bosons {h and H, which have 
scalar couplings to fermions at the tree level, and A with pseudoscalar tree-level couplings 
to fermions), all mix together in the presence of CP violation. The resulting three physical 
mass eigenstates have mixed CP parities. We denote these by hfi, 2,3 in order of increasing 
masses, i.e., Mhi < Mh 2 < Mhs- As a consequence of CP violation, all the three neutral 
Higgs bosons can now have tree-level couplings to pairs of and Z bosons. The Higgs 
couplings to the gauge bosons obey the sum rules 

3 3 

E 2 2 ^ ^ 2 2 / \ 

9HiW+w- = 9hw+w- 1 9HiZZ = 9hzz-i ( 1 - 1 ) 

1=1 i=l 

where Qhw+w- and Qhzz are the Standard Model couplings. Additionally, there is an 
important complementarity relation between the HiZZ and HiHjZ couplings: 


9HiHjZ — £ijk9HkZZ, 


( 1 . 2 ) 


where Eijk is the fully anti-symmetric Levi-Civita tensor. We have shown that the phe¬ 
nomenological consequences of (|l.lj) and (|1.2|) can be dramatic 0. In particular, (|1.1|) 
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implies that the HiZZ coupling can be significantly suppressed, thus raising the possibility 
that a relatively light Hi boson, with a mass even as low as 60 GeV, might have escaped 
detection at LEP 2. On the other hand, (|1.2|) leaves open the possibility of H 1 H 2 produc¬ 
tion at LEP 2, as discussed below. However, this brief discussion highlights the fact that 
a large part of the parameter space related to the Higgs sector must be re-explored in the 
presence of CP violation. 

In this paper, we investigate systematically the physics potential of the Standard 
Model Higgs boson searches at the Tevatron collider and the LHC for observing neutral 
Higgs bosons in the MSSM with explicit CP violation. We focus on Standard Model 
search channels, since the existence of a light Higgs boson with a significant coupling to 
W and Z bosons is a prediction of weak-scale supersymmetry. The present work goes 
beyond previous work on this subject m by incorporating the most important radiative 
corrections to the Higgs sector and performing a numerical analysis based on the most 
realistic simulations of Standard Model Higgs boson searches at the Tevatron and LHC. Our 
phenomenological study is performed mainly within the context of the CPX scenario 0, 
which was chosen to showcase the possible effect of CP violation in the Higgs sector of the 
MSSM. Given the limits from the LEP collider and existing simulations of the capabilities 
of the Tevatron and the LHC, we show that there are small regions of parameters in 
which all three neutral Higgs bosons escape detection. In these regions of parameters, the 
heavier Higgs bosons with the dominant coupling to the W and Z boson decay into Higgs 
boson with a tiny coupling to the W and Z boson, or into a lighter Higgs boson and a 
gauge boson. Outside these small regions of parameters, our analysis indicates, even in 
the presence of CP violation, the discovery of at least one neutral MSSM Higgs boson in 
a set of complementary detection channels. This provides a strong motivation to extend 
the experimental simulations to Higgs boson decay channels containing lighter Higgs boson 
states. For other choices of CP-violating parameters, we find that all three neutral Higgs 
bosons can have significant couplings to W and Z bosons, while being closely spaced in 
mass. This situation also requires careful investigation, since some signals may be affecting 
the background estimates for other signals. 

Our numerical analysis is performed using the code CPHDECAY |^, which is based on 
an extension of the code HDECAY Q to calculate the Higgs boson decay properties, and of 
the code cph to calculate the physical Higgs spectrum and mixing angles. The latter 

code includes the dominant one- and two-loop CP-violating contributions induced by the 
supersymmetric particles. At one loop, these contributions depend mostly on the phase of 
the trilinear mass parameters of the third-generation quarks, relative to that of the 
supersymmetric mass parameter fi*. At the two-loop level, a significant dependence on the 
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relative phase between /i* and the gluino mass parameter rrig also appears. Electroweak 
corrections are incorporated at the leading-logarithmic level. In a phase convention where 
the /i parameter is positive, the input parameters for the code are the top-quark mass 
the fi parameter, the soft supersymmetry-breaking masses Mq, Mf i, of the third-generation 
squarks, the third-generation soft trilinear couplings together with their respective 
phases, the gaugino masses and the phase of the gluino mass parameter rrig. Apart from 
its present application to hadron colliders, CPHDECAY may also serve as a useful tool for 
analogous studies at future linear colliders in or a muon collider [O, or for generic 


Higgs boson studies |T3 

Section ^ contains a complete discussion of the effective couplings of the neutral and 
charged Higgs bosons to quarks and of the Higgs self-couplings; full expressions are given 
in Appendix A. The one-loop corrected Higgs couplings are essential, as they determine the 
branching ratios of the Higgs bosons in the presence of CP violation. Section contains 
a discussion of the LEP results in the search for neutral Higgs bosons in the MSSM with 
explicit CP violation. Section ^ contains a numerical analysis of the different collider 
detection channels, showing graphically their complementary properties for the detection 
of Standard Model-like Higgs bosons in the MSSM with explicit CP violation. We restrict 
ourselves to the existing experimental simulations, identifying regions of parameters at 
small values of tan (3 and of the charged Higgs boson mass, for which none of the neutral 
Higgs bosons can be detected at the Tevatron or the LHC via the standard search channels. 
Our conclusions are presented in Section ^ Appendix A contains detailed formulae for 
couplings and Appendix B summarizes the information from detector simulations on which 
we base our analysis |2^,|4|,|2^,^ . 


2 Effective Higgs-Boson Couplings 

In order to study the collider phenomenology of the Higgs bosons in the MSSM with explicit 
CP violation, we must first calculate the couplings of the Higgs bosons to Standard Model 
and MSSM particles. In this section, we present explicitly those couplings which are the 
most important inputs to our numerical calculations. In particular, we review the effective 
couplings of the neutral Higgs bosons to quarks, and we present the effective Higgs self¬ 
couplings. We also give the corresponding effective couplings of the charged Higgs boson 
to the up and down quarks, even though they are not needed for the present analysis. 
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2.1 Effective Higgs-boson couplings to quarks 


Following our previous conventions throughout this paper [^], we define the physical scalar 
components of the Higgs superfields Hi and H 2 as $1 = iT 2 ^\ = (0?*, — 0r) $2 = 

(0^, 02); respectively, where T 2 is the usual Pauli matrix. The weak Higgs eigenstates are 
related to their mass eigenstates and through the orthogonal transformations: 


0? 


cos (3 — sin (3 = sin (3 + cos (3 , 

r’l + On Hi + i(^ cos (3 G^ — sin (3 O^i Hi ^ 

+ 02i Hi + i(^ sin /3 G® + cos f3 Osi Hi ) 


1 

P2 L 
1 

P2 


V2 


( 2 . 1 ) 

( 2 . 2 ) 


where G^ and G° are the would-be Goldstone bosons associated with the and Z bosons, 
respectively, and O is a 3-by-3 orthogonal matrix that describes the mixing of the neutral 
Higgs states in the presence of CP violation [^. 

Our starting point is the effective Lagrangian for the neutral Higgs-boson couplings 
to the u- and d-type quarks. We neglect quark-mixing effects, as their direct relevance 
to Higgs searches is not important .Q Moreover, we assume that there is negligible mixing 
between the different generations of squarks, and we include weak-interaction effects only 
at the leading-logarithmic level. 

Under these assumptions, the effective Lagrangian reads: 


_ po _ 


{hd + 5hd)ct)l* + 


bl* 


dn di 


^K4>\ + {K + 5hu)4>l urUl + h.c.. 


(2.3) 


where 5hd,u/hd,u and Ahd^u/hd,u are threshold radiative effects given by 
6hd 2a.« ^ . 


I'd 

hd 

^hd 

hd 

hu 

Shu 


dvr 


2q; 


* . * 


37r ^ 

2q; 




37r ^ 

2q; 




s * 

-m-y 

37r ^ 


2 2 1 |2\ 

) - 

15^2 \h\'^ , 

(2.4) 


^2^2 I l/ip), 

(2.5) 



(2.6) 

2 2 1 |2\ 


(2.7) 


We would like to emphasize that no other work includes the Ah;, corrections, and that 
these have important phenomenological implications, as we discuss below. 

^For a recent discussion of CP-violating Higgs-mediated flavor-changing neutral-current effects on low- 
energy observables, see Ref. []l^. 
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In (|2.4j) - (|2.7|) , ttg = g‘^/{4:7i) is the SU(3)c fine-structure constant, I{a,b,c) is the 
function 

nh]r\ {n / h\ hr ]n(h / r\ -\- nr]r\( r / n\ 

( 2 . 8 ) 


_ a 61 n(a/ 6 )-|- 6 cln( 6 /c)-|-acln(c/a) 

-/ ( Gj j (-) j C ) - 


(a — h){h — c){a — c) 
and 52 (with q = u,d) are the squark masses, with 
1 


cos2/3M| (2.9) 

-F(-) - M2 cos2/3M| (T| - 2QgSm'^ 6^) + 4m‘^ \Aq - Rqgi*\^ 

and QuiQd) = 2/3 (-1/3), T“ = -T^ = 1/2, Ru{Rd) = cot jS (tan/?), and sin^6'^ = 
1 — M^/M|. In ( ^73|) , the complex Yukawa couplings hd and hu are defined in a way 
such that the d- and n-quark masses are real and positive after the inclusion of radiative 
corrections, namely 

rrid 1 


hd = 

hu = 


vcos(3 1 + {6hd/hd) + (Ahrf/hrf) tan/? 
rriu 1 


( 2 . 10 ) 

( 2 . 11 ) 


nsin/? 1 -|- {6hu/hu) + (Ah„//iu) cot/? 

Substituting and (|2.1(]|) into (|2.3|) , the effective Lagrangian for the Hiqq couplings can 
be written in the general form: 

p Qw^q TJ - f 

^Hiqq — — ^ oTfY ^ ^ 


, 2M\^ ■ 1 

q=u,d 1=1 


[gkqq + ^9H,.qqh5] Q , 


( 2 . 12 ) 


with |§|,|5 


9%dd — 


Oii 


p 

9Hidd 


+ Im 

— Re 

— Im 


1 + Kd tan /? / cos /? 
Kd (tan^ /? -f 1) 


+ Re 


_ 1 + Kd tan /? . 


Ou , 


9%uu = Re 


tan (3 — Kd \ 

1 + Krftan/? J 

i^d \ _ 

1 + Kd tan /? / cos /? ’ 

1 \ O 


O^i + Im 

02i 


2 i 


P _ 
9HiUU 


Im 

Re 

Im 


1 + Ku cot /? / sin /? 
Ku (cot^ (3+1) 


+ Re 


. 1 + Ku cot (3 . 
cot (3 — Ku 


Osi , 


1 + Ku cot (3 


1 + Ku cot (3 J sin /? 


Osi Im 
Ou 


f^d ^ 

j 02i 

1 + Kd tan/? y 

' cos/? 

Kd tan/? / 

1 Ou 

1 + Kd tan (3 ) 

COS (3 


On 

1 + Ku cot (3 J 

sin (3 

Ku cot [3 ( 

1 02i 

1 + Ku cot (3 ) 

' sin;d 


(2.13) 


(2.14) 


(2.15) 


(2.16) 
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In the above, we have used the abbreviation: Hg = {Ahg/hq)/[l + {6hq/hq)], for q = u,d, 
and the fact that u- and d-quark masses are real and positive, so that Imm^ oc Im [hd + 
{5hd) + {Ahd) tan/3] = 0 and Imm„ oc Im [h^ + {Shu) + {Ahu) cot /3] = 0. 

Though our phenomenological analysis is restricted to the neutral Higgs-boson sector, 
we consider for completeness the charged-Higgs-quark sector. The effective Lagrangian 
describing the interactions of the charged-Higgs fields with u and d quarks is given by 


— — 


{hd + + Ahrf02 SrUl 

Ahuit>t + [K + 02 “R 0 R + h.c. 


(2.17) 


The quantities Shu^ and Ahu^ contain the respective threshold radiative effects that affect 
the charged-Higgs couplings to quarks. Though they are related to the 5hq and Ahq given 
above for the neutral Higgs sector, they differ by SU(2)^-breaking terms. 

To facilitate our presentation of the analytic forms of 5hq and Ahq, we hrst set up our 
conventions for the squark-mixing parameters and define some auxiliary functions. Thus, 
the left- and right-handed squark helds ql and qr are related to the physical helds qi ^2 
through the transformations: 


Ql = Cqqi + Sqq2 

Qr = -e^^'^Sqdi + e*'^’Cqg2 (2.18) 


where 6q = arg — Rqfi*) and Cq = cos 9q and Sq = sin6*g are squark-mixing angles. 
Analytic expressions for Cg and Sg may be found in Appendix A of p. We introduce the 
following short-hand form for the squark-loop integrals: 


I{i,j,m) = Imp), (2.19) 

where the indices i and j refer to the third-generation up- and down-squark mass eigen¬ 
states. 


With the aid of (|2.19|) , the following two auxiliary functions can be defined: 


Ii{m) 

l2{m) 


clsll{l,l,m) + slcll{2,2,m) + clcll{l,2,m) + slsll{2,l,m), (2.20) 
slcll{l,l,m) + clsll{2,2,m) -f slsll{1,2, m) + clcll{2,l,m). (2.21) 


Again, in the limit of neglecting mixing of generations, the threshold radiative effects related 
to the charged Higgs sector can be expressed in terms of (p.20|) and (|2.21|) as follows: 


Shd 

hd 


2as 

Stt 


m*-Adh{niq 


\hu? 

IOtt^ 


\h?h{h) , 


( 2 . 22 ) 
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^ (2.23) 

f: = (2-24) 

^ (2.25) 


Then, the effective Lagrangian for the charged Higgs-boson couplings to u- and d-type 
quarks is given by 


i^H±ud = ^{^^3H+udP- + T^d9H+udP+)d + h.c., (2.26) 

where = [1 =f 75 ]/ 2 , and 


9H+ud 


9H+ud 


cot^ (1 + Pu) - Ku 
1 + Ku cot P 

tan/3(l + p*) - R*^ 
1 + K*^ tan f3 


(2.27) 

(2.28) 


The quantities Ku,d are given after (|2.16|) , while Ku,d and pu,d in (|2.27|) and (|2.28|) are defined 
as follows: 

(Ahq/hq) _ {6hq/hq) - {5hq/hq) 


Kq = 


1 + {5hqlhq) ’ 
with q = u,d, respectively. 


PQ = 


1 + {Shq/hq) 


(2.29) 


It is apparent that the effective couplings of the neutral and charged Higgs bosons 
to quarks can be written in a form that depends entirely on the parameters Ku^d, Pu 4 and 
Pu^d- If is technically interesting to note that all these quantities are renormalization-scale 
independent up to the one-loop order in our resummation approach. More explicitly, Ku 4 
and Ru,d are proportional to the non-holomorphic radiative corrections {Ahu^d/hu,d) and 
{Ahu^d/hu,d), respectively, which are manifestly scale-independent. The quantities pu^d are 
measures of SU(2)£ breaking in the up- and down-Yukawa sectors. The parameters p^^d are 
also renormalization-scale independent, as a consequence of the SU(2)i gauge invariance 
of the original theory before spontaneous symmetry breaking. The scale independence of 
Pu,d can be also understood by simply noticing that the ultra-violet infinities, e.g., of 6hd 
and 6hd are equal, as they emanate from the charged and neutral Higgs components of 
the same gauge-invariant operator, which is QHiD in this case. The scale independence of 
these quantities has some analogies with that of Veltman’s p parameter that characterizes 
the isospin breaking in the SM gauge sector through the difference between the H^H^ and 
ZZ self-energies, that we do not develop further here. 
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Finally, we should bear in mind that the validity of the effective neutral and charged 
Higgs-boson couplings to quarks depend on two kinematic conditions: (i) the soft supersym¬ 
metry-breaking masses should be much larger than the electroweak scale and (ii) the ex¬ 
ternal momenta of the quarks and Higgs bosons, e.g., Pq^n, have to be sufficiently smaller 
than the soft supersymmetry-breaking mass scale Mgusv, so that they can be neglected 
when expanding the vertex functions in powers of p\h!■ 

2.2 Effective Higgs-boson self-couplings 

To determine the effective Higgs self-couplings, one needs to know the analytic forms of 
both the proper vertex and self-energy graphs. So far, there is complete information only 
for the latter contributions in the MSSM with explicit CP violation [^,^, while the former 
effects have been computed in |^, in an expansion of the effective potential up to operators 
of dimension 4. 

Here, we combine the above two pieces of information to obtain approximate analytic 
expressions for the trilinear and quadrilinear Higgs-boson self-couplings.^ Our analytic 
expressions are not limited to the MSSM case, but can be applied equally well to the 
general Two-Higgs-Doublet Model (2HDM) with explicit CP violation. 

To start with, we hrst write down the effective Lagrangian containing all operators 
of dimension 4: 


C = Ai(<f1$i)2 + A2(<h^<h2)' + A3(<hl<hi)(<F54>2) + A4(<hl<F2)(4>^<hi) 

+ A5(4>I<h2)' + A*($^<Fi)2 + A6($I<Fi)(<fI<F 2) + Aa<hl4>i) (<F^<Fi) 
+ A7(4>^4>2)(<hI<F2) + a; (<F^<F2) (4>^4>i), 


(2.30) 


where $^2 = (0p2) 0i2) two Higgs doublets, with their individual charged and 

neutral components defined in ( |2.1| ) and (|2.2|) . After (|2.1|) and ( p.2|) have been inserted into 
( p.30| ), the effective trilinear and quartic Higgs self-couplings may be cast into the form (in 
the unitary gauge): 



^<m<n=l,2,3 


^=1,2,3 


similar procedure has been followed in Ref. [^. 
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m<n<r<s=l,2,3 


AH 

mnrs i 


ij 


— On Omj 9l 

Km=l,2,3 


2HH+H- 

m 


and 

T^m+ 


— s^Ai + Cj^\2 + SpCpi^X^ + A 4 ) + 2s^c^ReA5 — 2s^c^ReA6 — 2s^c^ReA7. (2.34) 


The quantities gfjl, 9^^*^ , 9ijkii (p.33|) characterize the proper vertex correc¬ 

tions to the corresponding Higgs self-couplings. We observe that, exactly as is the case for 
in (|2.34|) , all these newly-introduced couplings depend on the radiatively-corrected 
quartic couplings Ai, A2, ..., Ay. Their analytic forms are given in Appendix A. 


■^4H+ 


3 Constraints on the CPX scenario from LEP searches 


The OPAL Collaboration |2^ has reported preliminary results on the search for Higgs 
bosons in the MSSM with explicit CP violation, taking the parameters to be those dehned 
in the CPX scenario IHI, 


AfsusY = Mq = Mt = Mb = 0.5 TeV, /i = 2 TeV, \At\ = \Ab\ = 1 TeV, 
\'^b\ = \'^w\ ~ 0.2 TeV, |mg| = 1 TeV, (3.35) 


where Msusy is the characteristic third generation squark mass scale, and are the 

bino and wino masses, respectively. 

OPAL presented exclusion regions in the M/y^-tan/? plane for different values of 
the CP-violating phases, assuming argA^ = argAb = argm^ = 0 cpx, with 0 cpx = 
90°, 60°, 30°, and 0°. Quite generally, there is no reason to expect the equality of the phases 
of the trilinear scalar couplings and the gaugino masses, nor to assume that 0 cpx < 90°, 
and our analysis departs from this assumption. 


To reproduce the results of the OPAL analysis, we rely on the combined ALEPH- 
DELPHI-L3-OPAL (ADLO) results for the ZLfgM(—> hh) channel [^| and the hA Ah 
channel [p^. Although the experimental hA analysis has been done for approximately 


equal values of the masses of the neutral Higgs bosons h and A, there is no major loss in 
efficiency when the splitting between these masses becomes larger |^, and therefore it 
is safe to apply those limits to a more generic set of masses Hi, Hj in the CP-violating 
scenario. 


Using the results of these experimental analyses, we have generated Figure 0, based 
on the CPX scenario with MgusY = 0.5 TeV. In this plot, the light grey area covers 
the theoretically allowed region of parameter space (consistent with electroweak symmetry 
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breaking), the medium grey region shows the exclusion from ZHi hnal states, the dark grey 
region is excluded by the search for Z* —>■ HiHj 46 hnal states, and the black region is 
excluded by both searches. 


When comparing the results in Fig. |l| with the OPAL results [^, we note that our 
limits on tan/3 are somewhat stronger than the OPAL limits. This can be attributed 
to the fact that our analysis is based on combined ADLO limits obtained from all four 
experiments The results of our analysis show also a somewhat better reach than the 
OPAL one in the HiHj production channel. This can again be attributed to our use of the 
combined results of all four experiments and the fact that, following the above discussion, 
we have assumed no deterioration of the H^Hj signal for ^ We expect that our 

CPX estimates will be quite similar to the hnal LEP combined results. 

We see from the cases displayed in Fig. that the case of vanishing phases is most 
severely constrained by the LEP data. Most of the coverage arises from the ZHi or H 1 H 2 
processes. The appearance of two ‘hngers’ in the ZHi coverage at large tan/3 for vanishing 
phases arises from the shape of the LEP exclusion curve, which is distorted by a marginal 
excess in the region rrih ~ 80 — 90 GeV [^]. The case of phase 30° is very similar, and the 
60° case has analogous features. 


For signihcant values of the phases, there are regions of parameters in which the 
lightest neutral Higgs boson is very weakly coupled to the Z gauge boson, and light enough 
for the heavier neutral Higgs states to decay into a pair of Hi bosons. In these regions 
of parameters, the heavier Higgs states decay into a hnal state containing four b quarks. 
Therefore, the dominant production and decay modes contain 6 jets in the hnal state. 
The current experimental strategy is to employ the standard 4-6 analysis by forcing (with 
suitable jet dehnitions) the 6-jet topologies into 4-jet ones. This leads to a very low efficiency 
for the real 6-jet signal, and hence a signihcant region of CPX parameter space where light 
Higgs bosons may exist but are not excluded by the current analyses. In our analysis 
of the MSSM parameter space excluded by Higgs boson searches at LEP, we have not 
attempted to simulate the signatures associated with the decay of heavy neutral Higgs 
bosons into lighter ones. Nonetheless, we are hopeful that a dedicated 6-jet analysis could 
cover this region. Returning to Fig. ^ the upper two panels (90°, 60°) demonstrate that 
LEP cannot exclude the presence of a light Higgs boson at tan/3 ~ 3 — 5, Mh^ < 60 GeV 
and tan /3 ~ 2 — 3, Mh^ ~ 40 GeV, respectively, in good qualitative agreement with the 
OPAL results. Using the ADLO data in the Zifj(—> 66) channel, a Higgs boson with the 
Standard Model coupling can be excluded if <115 GeV, and even a low-mass Higgs 
boson can be excluded if \gHiZz\ ~ -22 relative to the Standard Model coupling. 

In Figs. |]and ^ we have presented a generalization of the above analysis to unequal 
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LEP Exclusion for CPXq ^ 
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Mhi (GeV) 

Figure 1: Approximate LEP exclusion limits in the plane for various CPX 

scenarios, using combined LEP results. The light grey covers all the region of parameter 
space that is consistent with electroweak symmetry breaking, the medium grey shows the 
exclusion from e+e“ — ZHi, the dark grey shows the region excluded by Z* HiHj 46 
searches, and the black region is excluded by both searches. 
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LEP Exclusion for CPX^ ^ 
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Mj^+ (GeV) 

Figure 2: LEP exclusion limits in the CPX benchmark scenario in the Mh+Xqh j3 plane 
for different values of the phases arg{A) and arg{mf) of the trilinear couplings At^b and 
the gluino mass parameter, respectively. The four panels show the results for (0,0); (90,0); 
(90,90) and (135,90) degrees, respectively. The light grey region is disallowed theoretically, 
the medium grey region is excluded by the absence of ZHi, -hatched region by the absence 
of H 1 H 2 , 135°-hatched region by the absence of ZH 2 , and the horizontally-hatched region 
by the absence of HiH^. 
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Figure 3: Panels (a) and (b) are similar to Fig. but for CP-violating phases (60,60) 
and ( 140 , 140 ) degrees. Panels (c) and (d) show two special scenarios 3HGC1 and 3HGC2, 
which are particularly challenging for all colliders: 3HGG1: Msusy = 0.5 TeV, A = 1.95, 
p. = 2.4, arg{A) = arg(mg) = 75°; 3HGG2: Msusy = 0.5 TeV, A = 2, p = 2, arg{A) = 
arg(mg) = 145°, M~ = 0.5 TeV. The shadings and hatchings have the same significances 
as in Fig. 
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values of the CP-violating phases of the trilinear mass parameters f, and the gluino mass. 
The exclusion plots are presented in the Mji/+-tan/3 plane. For the case of equal values 
of the phases, these hgures present analogous information to that presented in Fig. 
Note the appearance of hngers for m~) = (0°, 0°), as well as uncovered regions near 

tan/3 ~ 4 — 5 and Mjj± ~ 125 — 140 GeV for (90°, 90°), and near tan/3 ~ 2 — 3 and 
Mh± ~ 105 — 130 GeV for (60°, 60°). Moreover, Figs. | and ^ show distinctively the 
regions covered by the different channels studied at LEP. The light grey regions of these 
hgures are excluded theoretically, the medium grey regions are excluded by the absence of 
ZHi, the regions hatched at 45° by the absence of H 1 H 2 , the regions hatched at 135° by 
the absence of ZH 2 , and the horizontally-hatched regions by the absence of HiH^. 

Fig. H is a study of the interplay of phases for the trilinear parameters At^i, and 
the gluino mass parameter. The appearance of substantial phases to the parameters can 
signihcantly modify the collider phenomenology (the phase of At providing the dominant 
effect). The addition of a phase to the gluino mass parameter, however contributes to GP- 
violating effects only at the two-loop level and its effect is comparatively smaller, as can be 
seen in Fig. H. For the case of axg{At^b) = 90°, it is clear from Fig. that, for tan/3 ~ 4-5, as 
the Higgs Hi becomes lighter, it starts decoupling from the Z. There is a signihcant region 
of parameter space where, although heavier, H 2 is in the kinematic region accessible to 
LEP and couples relevantly to the Z gauge boson. For Mh+ ~ 130 GeV, Hi becomes light 
enough that the decay H 2 —> HiHi dominates over H 2 —>■ bb, and therefore H 2 detection 
becomes difficult. Finally, for arg(Ht f,) = 135° and arg(mg) = 90°, and for moderate or 
large values of the charged Higgs mass, all neutral Higgs bosons rapidly become sufficiently 
heavy to be out of the reach of LEP for moderate and large values of tan /3. 

In Fig. ^ we display results for the GPX scenario, for equal values of the At^b and 
gluino mass phases equal to 60° and 140°. The region of parameters left uncovered by LEP 
for low values of tan/3 and of the charged Higgs boson mass, for phases equal to 60°, is 
seen more clearly than in Fig. It is also apparent that, as happens for a GP-violating 
phase equal to 90°, close to the region left uncovered by LEP, there is a large region of 
parameters covered by the ZH 2 channel. 

For phases equal to 140°, we hnd a peculiar behavior of the covered region of param¬ 
eters for large values of the charged Higgs mass, which can be traced to the behavior of the 
Hi mass shown in Fig. |^. For small and moderate values of tan/3, it has the standard be¬ 
havior of the MSSM without explicit GP violation, increasing with tan /3. For tan /3 above 
about 10, however, it starts decreasing due to the effect of radiative corrections from the 
sbottom sector, which involve the bottom Yukawa coupling. The bottom quark Yukawa 
corrections are screened for sufficiently large values of tan /3, causing the mass of the Hi to 
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Figure 4: The non-monotonic behavior of Mh^, the mass of the lightest neutral Higgs boson, 
in the CPX scenario with arg(At) = arg(m~) = 14(T for different values of Mh+. The dip 
for 15 < tan < 30 is due to the behavior of the radiative corrections from the sbottom 
sector, which are screened at large tan f3, as discussed in the text. 


increase again and eventually become larger than the LEP kinematic limit. 

This screening phenomena occurs because the sbottom-induced radiative corrections 
are complex. At large values of tan/?, the bottom-quark Yukawa coupling h^ is approxi¬ 
mately related to the bottom mass and to the Yukawa corrections Xh^ by 


^ ruf, tan/5 

^ V 1-|-(Ahfe/hfe) tan/5 


The negative corrections to the Higgs mass squared coming from sbottom loops are pro¬ 
portional to the fourth power of the modulus of hh- For equal values of the phases of the 
gluino and of the trilinear mass parameter At, i.e. arg(At) = axg{mf) = 0cpx, the modulus 
of hb is given by 


mb tan (3 

\ll + 2\Xhb/hb\ cos(0cpx) tan/5 -1- (|Ahfe//ib| tan/5)^ 


(3.37) 


For cos(0cpx) > 0, the above expression is a monotonically increasing function of tan/5, 
meaning that the negative corrections induced by sbottom loops become larger for larger 
values of tan /5. 

On the other hand, for negative values of cos(0cpx), the expression for \hb\ has a 
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maximum at a value of 


tan/3”“ = 


(3.38) 


|A/ib//ib| cos(0cpx) 

meaning that the negative effects on the Higgs mass induced by sbottom loops are most 
pronounced at intermediate values of tan/3, and are reduced again at large tan/3. For 
the particular case under consideration, we have |Ahb|/hfc ~ 1/20 and 0cpx = 140°, with 
cos(0cpx) — l/\/2, while \hh\ ~ 1 stays well within the perturbative range. Thus, lower 
values of the Higgs mass are obtained for tan /3 ~ 25 than for either smaller or larger values 
of tan/3. This feature is seen clearly in Fig. and explains the stronger coverage of LEP 
at this intermediate region of tan/3, as seen in panel (b) of Fig. 0. 


Apart from the above property, the only salient feature is that in the region between 
the one covered by LEP in the ZHi channel and that covered by LEP in the H 1 H 2 channel. 
Hi becomes light but couples very weakly to the Z boson. 


The last two scenarios considered in Fig. are interesting since, for moderate values 
of the charged Higgs mass of about 150 GeV, the three neutral Higgs bosons become quite 
similar in mass and share approximately equal couplings to the weak gauge bosons. We call 
these tri-Higgs-gauge coupling (3HGC) scenarios, which are denoted in short as 3HGG1 
and 3HGG2. As we discuss below, the searches at hadron colliders become particularly 
difficult in 3HGG1 and 3HGG2. The excluded region for moderate and large values of the 
charged Higgs-boson mass and large values of tan/3 in the scenario 3HGG2 is related to 
the Higgs-sbottom coupling effect discussed above. 


To summarize the results from LEP searches, the phenomenology is mainly sensitive 
to the values of gUiZZ and the kinematic limit on the Higgs boson masses. The exceptions 
are the cases when the decay H 2 —> HiHi becomes relevant at small tan/3 and when 
radiative corrections extend the LEP coverage for certain choices of phases </>cpx at large 
tan/3. A Higgs boson with the Standard Model coupling can be excluded if M//. < 115 
GeV, and even a low-mass Higgs boson can be excluded if \gHiZz\ ~ -22 relative to the 
Standard Model coupling. Note that the application of the ZHi and HiHj results to the 
GPX scenarios treats each potential signature independently. Any possible complications 
from one potential signal effecting the background estimate for another (and vice versa) 
are ignored. We are not aware of any experimental analysis of such a situation. 


4 Higgs Boson Searches at Hadron Colliders 

Experiments at the Tevatron and LHG will probe the MSSM Higgs sector beyond the 
kinematic and dynamical reach of LEP. At these hadron colliders, we have considered 
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several channels for Higgs boson Hi searches, where Hi stands for -Hi, 2 , 3 : 

(a) tiHi[-^ bb) at the LHC 

(b) W/ZHi{^ bb) at the Tevatron 

(c) WW Hi{^ T+T-) at the LHC 

(d) gg ^ Hi ^ 77 at the LHC. 

In the absence of reliable experimental simnlations, as in the case of LEP, we have not 
modihed onr analysis to acconnt for the decay Hj HiHi —>• 46 (j > i), shonld it 
become relevant. The search channels (a) — (d) are considered to be the most promising 
for observing a Standard Model Higgs boson, and have been the most thoronghly analyzed 
and simnlated. Since the MSSM Higgs sector has a deconpling limit, where there is a 
light Higgs boson with properties almost indistingnishable from a SM Higgs boson, it is 
appropriate to apply these SM analyses to the MSSM case. However, one wonld have to 
be carefnl in interpreting a signal in some of these channels. For example, signatnre (d) 
may arise when Hi is replaced by a psendo-Goldstone boson, as in Technicolor models. 
Signatnre (a) may arise, for example, when Hi is replaced by a top-pion, as in a Topcolor 
model. Only observation of channels ( 6 ) and (c) wonld clearly indicate that the scalar Hi is 
associated with EWSB. Signatnre (a) may indicate that the scalar responsible for EWSB 
also generates fermion masses. Nonetheless, if signatnre (a) or (d) were observed with 
the rate expected in a weakly-conpled theory, it conld be considered strong evidence for a 
fnndamental Higgs scalar or scalars. In the present work, we do not analyze those signatnres 
associated with non-SM-like Higgs bosons, snch as charged Higgs or psendoscalar Higgs 
bosons, since their interpretation wonld be ambignons. 

Fignre shows the coverage of the M^/^-tan/? plane by the Tevatron W/ZHi{—>- bb) 
search (45° lines) and the combined LHC coverage for the tiHi{-^ bb), WW —> Hi{-^ t^t~) 
and gg Hi{—>- 77 ) searches (135° lines) for the same valnes of the CP-violating phases as 
chosen in Fig. |I[ For the Tevatron, we show the 3a evidence coverage with 5 fb”^, while, 
for the LHC, we show the 5a discovery coverage for the 77 and bb channels with 100 fb“^ 
and the channel with 30 fb“^. The previons LEP 95% C.L. exclnsion is also inclnded 
(medinm grey) snperimposed on the theoretically allowed region (light grey). While the 
Tevatron and LHC searches are adeqnate for extending the coverage into the large-tan/3 
region, the region of small Mjj^ and low tan/3 for 0 cpx = 90°, 60° remains nncovered. 
Althongh the persistence of these regions can be clearly traced to the decay of the heavier 
Higgs bosons into lighter ones, one conld inqnire why the LHC cannot see the lighter Higgs 
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boson in the tiHi and/or the glnon fusion channels. The reason is that, in the same region 
of parameters, Hi couples weakly to the top quark, as well as negligibly to W and Z bosons. 
Therefore, not only is the ttHi production channel suppressed but the loop-induced gluon 
fusion production and decay into photons is suppressed as well. Furthermore, since tan (3 
is small. Hi would likely not be observed in the pseudoscalar Higgs channel bhHi{-^ bb). 


If, contrary to our expectations, a dedicated LEP analysis cannot cover this prob¬ 
lematic region, there still remains the possibility of observing H 2 —HiHi decays (when 
H 2 is otherwise SM-like) at hadron colliders. Since H 2 is SM-like, there still is a substan¬ 
tial production cross section. Dedicated searches should be able to identify 46 hnal states 
with high efficiency and substantially less background than for the 26 case. Various cases 
of Higgs pair production, such as gg —>■ H —y hh/AA, have been considered in previous 
studies for the Tevatron and LHC, but mostly at the parton level |^,|^. The studies 
would have to be updated for the case of explicit CP violation. Searches for W/ZH 2 {—y 46) 
or ttH 2 {-^ 46) may be more promising, if only because of the reduced backgrounds from 
requiring the leptonic decay oi a W boson. 


Outside these small regions of parameters, the search for neutral Higgs bosons at the 
Tevatron and LHC colliders are complementary and sufficient to cover the full parameter 
space in one or several channels. The Tevatron does not provide any additional coverage 
alone over all three LHC search channels. However, if we just focus on the W/ZHi and 
WW — > Hi channels - the only two which are dependent on the coupling gntww ^ then 
both are nearly sufficient to cover the entire parameter space: for 90°, Mh^ ~ 80 GeV and 
tan/3 ~ 10, a small region is left uncovered without tiHi. This complementarity reflects 
the different production and decay channels that are being used in the search for neutral 
Higgs bosons, and is similar to that in the CP-conserving scenario ||^. For instance, there 
are regions of parameter space where either the 66 or the decay branching ratio of 

the lightest neutral Higgs boson is suppressed with respect to the Standard Model ones, 
while the other neutral Higgs bosons are heavy and weakly coupled to the gauge bosons. 
These suppressions tend to occur for large values of tan The large values of /i chosen 
in our examples maximize the CP-violating effect, but also determine that for large values 
of tan/3 the quantum corrections to the bottom Yukawa coupling are large and sizeable, 
causing a visible displacement of the regions of suppressions of each coupling. Since, for a 
light enough Higgs, these two branching ratios tend to be very important, a suppression of 
one of them tends to enhance the other. We discuss specihc examples below. For smaller 
values of /i, both the 66 and couplings can be suppressed simultaneously and then the 
77 decay branching ratio of the lightest Higgs boson tends to be enhanced. An analysis of 
these cases can be found in Ref. . 
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Figure 5: Approximate Tevatron/LHC discovery and LEP exclusion limits in the Mh^- 
tau/S plane for the CPX scenario with both phases set to: (a) 90°, (h) 60°, (c) 30°, and 
(d) 0°. The reach of the Tevatron W/ZHi{-^ bb) search is shown as 45° lines and that of 
the combined LHC search channels as 135° lines. The combined LEP exclusion is shown in 
medium gray, superimposed on the theoretically allowed region in light grey. 
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4.1 Detailed Analysis of the different scenarios 


In Figs. HI, we present a detailed analysis of the regions covered by the relevant search 
channels at the Tevatron and at the LHC. In each of these Figures, the four panels display 
the regions for which Higgs bosons may be discovered at the 5-a level at the LHC via 
the tiHi channel with 100 fb“^ (two darkest shades of grey) and 30 fb“^ (135° lines) 
the weak-boson-fusion channel using 30 fb“^ (two darkest shades of grey) and the gluon 
fusion channel with Hi —>• 77 using 100 fb“^ (two darkest shades of grey) and 30 fb“^ 
(135° lines), as well as the region for which the Tevatron can provide 3 - 0 " evidence with 
5 fb“^ (two darkest shades of grey) and 2 fb“^ (135° lines) of integrated luminosity. In 
each case, the darker shade of grey extends beyond the LEP exclusion. The 3 -ct coverage 
with 5 fb“^ at the Tevatron is very similar to the region that would be obtained for 5-a 
discovery with 15 fb“^, or for 95% C.L. exclusion with 2 fb“^. The reason why we have only 
analyzed 30 fb“^ in the case of weak-boson fusion is that a realistic analysis exists only for 
moderate instantaneous luminosity, and detection in this channel might be compromised 
in the different environment created by high-luminosity LHC running. For the Tevatron, 
we chose to present results that can simultaneously indicate the Higgs discovery potential 
under a conservative assumption of the total integrated luminosity (5 fb“^), for which 
only a 3-cr evidence may be obtained in the regions not excluded by LEP. This can be 
directly compared to a more optimistic assumption (15 fb“^), for which discovery may be 
possible well beyond the region of parameters covered by LEP. Finally, in all the scenarios 
studied, there are regions where at least two Higgs bosons are close in mass and have 
similar strength couplings to W/Z bosons. Our prescription is to add the signals from 
Hi —>■ bb and Hi with no degradation if the mass difference is less than 5 GeV 

(which is substantially smaller than the expected mass resolution), otherwise the signals 
are treated separately. We discuss this more later. Because of the precision electromagnetic 
calorimetry expected at the LHC, we do not combine Hi 77 signatures, though this is 
of little practical consequence. 


Figure ^ demonstrates the case of vanishing phases. Because the CPX scenario has 
large values for ji and m^, there is a large value for the supersymmetric loop corrections 
to the &-quark Yukawa coupling from (|2.4| , p2.5| ). Furthermore, since sgn{fim~) > 0, the b 
quark Yukawa coupling is decreased relative to the r lepton Yukawa coupling. The right 
lobe of the tiHi coverage (a) arises from Hi = Hi, with a transition to H 2 and then H^ 
in the left lobe. The suppression of the bb coupling relative to , in conjunction with 
the sharing of the gHat coupling between (in this case) two Higgs bosons, results in the dip 
near Mu± ~ 150 GeV. For this case without explicit CP violation, that region corresponds 
exactly to ~ where is the maximal value of the Standard Model-like Higgs 
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boson in the deconpling limit. Because of the enhanced coupling relative to bb, the 

WW —>■ h channel (c) does not exhibit this phenomena. For tan/3 > a few, the W/ZHi 
channel (b) has a similar behavior as tiHi. The limited coverage from the h —> 77 channel 
(d) is a common feature of all the scenarios we have studied. Even when gUiWW ~ 1, 
a suppression of BR(h — 77 ) occurs because is enhanced over the SM value. The 
large values of |/r| and 1^41 necessary to cause substantial effects typically increased the off- 
diagonal elements of the 3x3 Higgs squared-mass matrix relative to the diagonal elements. 
Note that the WW —^ channel with 30 fb“^ is sufficient alone to cover the 

M^i-tan/? plane. With the same integrated luminosity, neither the tiHi nor the Hi —> 77 
channels even extend beyond the LEP coverage. Of course, the former is a discovery region, 
whereas the latter is only 95% C.L. exclusion. 

Another possible variant on the CPX scenario is to introduce CP violation to the 
Higgs sector through a substantial phase for m~, (jy^ = 90° (not shown). The primary effect 
is to decrease the suppression of gHiU relative to gHirr, but, as discussed before, the gluino 
phase effects are quite mild. 

The introduction of a phase for At = At, allows for CP-violation effects at small 
values of tan/3, as shown in Fig. ^ As discussed previously, it is now possible for the 
lightest Higgs boson Hi to become quite light and have a negligible coupling to vector 
bosons. Furthermore, the decay H 2 HiHi becomes kinematically possible, introducing a 
region that remains uncovered by all experiments if one considers only the detailed analyses 
performed so far. Another region of difficult coverage appears in the tiHi{-^ bb) channel, 
for a charged-Higgs mass close to 125 GeV. This did not occur for vanishing phases, and 
only arises now because of the transition away from the region where H 2 —> HiHi occurs 
when H 2 has a substantial coupling to the Z boson. This does not occur for weak boson 
fusion in the Hi tt channel because of the enhancement of the r coupling over that 
of the b. In this region of parameters, the coupling of Hi and H 2 to the top quark and 
weak bosons is somewhat suppressed, though H 2 is relatively more strongly coupled to 
these particles than Hi. Due to hnite quantum corrections, the ratio of branching ratios 
BR(ifi —> 66 )/BR(ifi ^ tt) is suppressed with respect to the Standard Model value. The 
substantial increase in BR(if 2 —^ tt) means that H 2 —>■ tt can still be observed, even 
though its vector-boson coupling is suppressed. Indeed, H 2 is also visible in the channel 
tiH 2 in a small region of parameters, for somewhat larger values of the charged Higgs mass. 

We show in Fig. || the effect of including a gluino phase as well, which, as mentioned 
before, does not change signihcantly the picture, apart from the fact that small regions of 
uncovered parameter space appear. Those uncovered regions appear for a charged-Higgs 
mass around 130 GeV (160 GeV) because two of the Higgs bosons, H^, and H 2 {Hi and 
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Figure 6: Coverage of the plane in the CPX scenario with Msusy = 0.5 TeV 

and for the phases shown at top. The light grey region is theoretically excluded, and the 
two medium grey regions are excluded by LEP . The two darkest shades of grey are (a) the 
5-a discovery region at the LHC using ttHi bb (100 fb~^); (b) the 3-a evidence region at 
the Tevatron using WfZHi bb (5 fb~^); (c) the 5-a discovery region at the LHC using 
WW Hi ^ T^T~ (30 fb~^); and (d) the 5-a discovery region at the LHC using Hi 77 
with 100 fb~^ of luminosity. The 135° lines in panels (a), (b), and (d) show the coverage 
for 30, 2, and 30 fb~^ of data, respectively. This baseline model shows the behavior when 
no CP violation is present, i.e., the phases of At = At, and m~ are both zero. 
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H 2 ) are separated in mass (greater than 5 GeV apart) and share equal couplings, of order 
a half of the SM values, to the weak bosons and the top quark. 

An interesting feature appears for the case displayed in Fig. |^. For moderate and large 
values of the charged-Higgs mass and large values of tan/? ~ 30, the lightest Higgs-boson 
mass is about 123 GeV. For these values of the Higgs mass, discovery in the weak-boson 
fusion process demands a slightly larger branching ratio into r leptons than in the SM. 
This occurs for smaller values of the charged-Higgs mass, but as the charged-Higgs mass 
increases the branching ratio into r leptons becomes smaller than the SM one, and remains 
like this until the decoupling limit. The bb decay branching ratio has a different behavior, 
as we see in the panel showing the expected Tevatron coverage. The mass Mh^ is about 
123 GeV in the large-large-tan/? limit, but ^ bb) is slightly larger than in the 

SM, explaining the appearance of coverage in that region; the lack of coverage starting near 
Mh± 235 GeV and tan /? ~ 20 reflects the difficulty of establishing evidence for a SM 
Higgs boson of mass 123 GeV. As usual, for larger values of the charged-Higgs mass, the 
77 mode becomes relevant as a discovery channel, since the search for a Standard Model 
Higgs boson is efficient. Note also the disappearance of an uncovered region associated 
with H 2 —> HiHi. 

The case of phases equal to 60°, displayed in Fig. |^, is similar to the case of 90°, 
although the asymptotic value of the Hi mass is somewhat lower. Most of the features of 
the Higgs searches in these two cases have similar behaviors, though occurring for different 
parameters, so we do not expand on them further. 

Our hnal example for the GPX scenario has phases equal to 140°. In this case, 
unlike the other examples, there is a relative suppression of the r branching ratio of Hi 
in large fractions of the parameter space, which explains the large regions of parameter 
space uncovered by the weak-boson-fusion process. This search channel becomes efficient 
whenever the Higgs mass is at the limit of the LEP reach, since in this case only a lower 
decay branching ratio is sufficient, as seen in Table 3 of Appendix B. For slightly shifted 
values of Mh± and tan/3, instead the bb coupling is suppressed relative to The 

complementarity of these channels is clearly illustrated in panels (b) and (d). 

One interesting phenomena that arises in the GPX scenario, which is experimentally 
challenging, is the appearance of regions uncovered because of H 2 HiHi decay. As 
noted previously, there are regions where two Higgs bosons share the coupling io W/Z 
bosons and ti, but are separated in mass. In this ‘transition region’ from one Higgs boson 
being most SM-Iike to another one, the effective signal is halved. In the MSSM without 
explicit GP violation, these Higgs bosons are nearly degenerate in mass, and they are most 
likely indistinguishable from two Higgs bosons with the experimental resolution. After 


24 


including explicit CP violation, the mass splittings can be significantly larger, and all three 
Higgs bosons can share the WjZ coupling. We present here only two example scenarios that 
exhibit this phenomena - many more are possible. Under our prescription for adding signals 
occurring at masses within 5 GeV, the last two scenarios analyzed in Figs. |T^ and show 
small regions of parameters uncovered by any of the colliders. The size of these uncovered 
regions may be underestimated. Since there are potentially 3 Higgs bosons yielding a 
similar signature in these regions, and one may have to rely on experiment to normalize 
the background, the Higgs signal may be entirely washed out. Consider, for example, the 
expected number of signal and background events for the r+r” channel shown in Table 3. 
Since SjB ^ 1, the accumulation of signals in nearby bins could substantially increase the 
background estimate. Clearly, the possibility that Higgs bosons may provide backgrounds 
for the other Higgs bosons should be understood and analyzed in more detail. 


On the other hand, one consequence of sharing the coupling to W and Z bosons is that 
all three neutral Higgs bosons may be observable at hadron colliders in the WjZHi and/or 
WW —> Hi channels. This is an exciting possibility, since it would indicate that a 2HDM 
without explicit CP violation is inadequate to describe the data. We find no such overlap 
of 5-(T signatures from all three neutral Higgs bosons in these channels for the scenarios 
presented here. Furthermore, using the expected reach of the ATLAS collaboration for 
the gg —> H^m ZZ* ^ M channel [^, we have checked that there is almost no coverage 
beyond LEP for observing even a single Higgs boson. 


5 Conclusions 


We have presented in this article a phenomenological analysis of Higgs boson searches in 
Standard Model channels at the Tevatron collider and the LHC, for the case of the MSSM 
with explicit CP violation. We have also provided analytical expressions for the effective 
couplings of the neutral and charged Higgs bosons to fermions and to gauge bosons and the 
Higgs-boson self couplings within the MSSM with explicit CP violation. These expressions 
have been incorporated in the program CPHDECAY , used to perform our analysis. 


After considering the LEP limits and comparing our results with the published LEP 
analyses, we have analyzed the reach of the Tevatron collider and the LHC in the Standard 
Model search channels for a neutral Higgs boson in the CPX scenario, as well as in other 
interesting scenarios. Our study was motivated by the fact that the latest LEP data prove 
insufficient to exclude a light MSSM Higgs boson in the CPX scenario, with a mass smaller 
than 60-70 GeV. Further coverage of the MSSM parameter space with CP-violating phases 
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Figure 7: Same as Fig. but for CP-violating phases (arg(At^b), arg(mg)) = 
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Figure 8: Same as Fig. but for CP-violating phases (arg(At^b), arg(mg)) = (00°,90°j. 
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Figure 9: Same as Fig. but for CP-violating phases (arg(At^b), arg(mg)) = (^135°,90°j. 
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Figure 10: Same as Fig. but for CP-violating phases (arg(At^b), arg(mg)) = fB0°,60°j. 
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Figure 11: Same as Fig. but for CP-violating phases (arg(At^b), o,rg(mg)) = (^140°, 140°^. 
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Figure 12: Same as Fig. but for the 3HGC1 scenario defined in Fig. |^. Here, a region 
remains uncovered due to dilution of the gUiVV (V = W, Z) coupling among three Higgs 
bosons. 


31 












CO. 

c 

cd 


3HGC2q5 145® 145® 

100 125 150 175 200 225 250 275 100 125 150 175 200 225 250 275 


40 

30 


20 - 


10 - 


5 - 


2 

40 

30 


20 - 


10 - 


5 - 



p-I-TTT 


I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



J 


WW^R 






1111 1111 1111 1111 1111 1111 1111 1111 1111 1111 


. 


I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I H 





100 125 150 175 200 225 250 275 100 125 150 175 200 225 250 275 


40 

30 

20 

10 

5 


2 

40 

30 

20 

10 

5 


2 


Mj^+ (GeV) 

Figure 13: Same as Fig. \^, but for the 3HGC2 scenario defined in Fig. |^. ^45 in Fig. \l^, a 
region remains uncovered due to dilution of the gUiVV coupling among three Higgs bosons. 
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might still be possible with LEP data if searches for all the possible decay channels of the 
Higgs bosons were optimized. Our analysis for the Tevatron and the LHC is based on 
existing experimental simulations of the search for neutral Higgs bosons at these colliders, 
and could also be extended if more simulations were available. 

We have shown that, in general, there is complementarity between the various searches 
at the two colliders for neutral Higgs bosons with Standard Model-like properties. However, 
there are still small regions of parameter space, for small values of the charged Higgs boson 
mass and moderate values of tan (3, in which none of the three neutral Higgs bosons can 
be detected with a high statistical signihcance. In these regions of parameters, one of two 
phenomena occurs. First, the neutral Higgs boson with dominant couplings to the W and 
Z bosons can decay predominantly into channels which contain either two neutral Higgs 
bosons, or a neutral Higgs boson and a Z boson. The lighter Higgs boson has only feeble 
couplings to the W and Z bosons and top quarks, and escapes detection both at LEP 
and the hadron colliders. Secondly, all three neutral Higgs bosons can share the coupling 
to W and Z bosons and the top quark, resulting in three marginal signal excesses. The 
masses of the Higgs bosons may be uncomfortably close, leading to the possibility that one 
“signal” is a background to another. Detailed experimental simulations of these situations 
should be performed in order to decide the detectability of Higgs bosons in these regions 
of parameters. 

Apart from these regions, the Standard Model search channels at the LHC are ade¬ 
quate for discovering a neutral Higgs boson over the range of supersymmetric parameters 
considered here. The bb) and WW Hi{-^ t^t~) channels are most important. 

The gg ^ Hi ^ 77 channel is greatly suppressed, and provides no additional coverage. 
This is a generic feature of the CPX and related scenarios where the 6 -quark Yukawa cou¬ 
pling can be enhanced over its Standard Model value. However, the Tevatron W/ZHi{-^ bb) 
and the LHC WW Hi{-^ t^t~) channels fall in a special category for two reasons. First, 
the production of the Higgs boson through a tree-level coupling to the W and Z boson 
coupling demonstrates that this Higgs boson contributed to the Higgs mechanism. Sec¬ 
ondly, radiative corrections can enhance or suppress the 6 -quark Yukawa coupling, while 
not affecting the r-lepton Yukawa coupling, so that BR{Hi bb) and BR(iLj — t~^t~) 
can be complementary to each other. We have presented results for 3-a evidence with 5 
fb“^ of luminosity at the Tevatron, and 5-a discovery with 30 fb“^ of luminosity for the 
T~^T~ signature at the LHC and with 100 fb“^ for the bb and 77 channels. Our results 
indicate that the W/ZHi{-^ bb) and WW Hi{-^ t^t~) channels alone do not cover the 
entire range of supersymmetric parameters, even excluding the difficult regions highlighted 
previously. 
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The restriction of WW Hi{-^ t^t~) results to only 30 fb“^, which corresponds to 
accumulating data for a few years at a moderate instantaneous luminosity, is motivated by 
the potential deterioration of detection efficiency for this signal at a higher instantaneous 
luminosity. It is notable that, when comparing the reach with only 30 fb“^ of data for 
all search channels, the bb) channel provides very little coverage beyond LEP - 

though the LHC coverage is for discovery, whereas the LEP coverage is for exclusion - and 
the gg Hi{-^ 77 ) channel provides no discovery potential. It is clear that the vector- 
boson fusion channel should be studied more thoroughly with detailed detector simulations, 
including the potential complications raised here regarding several marginal signals in a 
similar kinematic region. 
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A Proper Higgs-boson self-couplings 


We exhibit here analytic expressions for the proper Higgs-boson self-conplings gfjl, , 

dfjkiJ ■ The Latin indices = 1,2,3 attached to the proper Higgs self- 

conplings correspond to weak states 0i, (j )2 and a. 

First, we list the proper trilinear self-conplings gf^^^ and gf^-. 


^ — 2 s^C/3Ai -|- c^As — s'^cgX^ — 2s^c^ReA5 -1- s_(3(s^ — 2c^) ReAg -f- spc^^KeXj, 

gHH H _ 2 s^c^A2 -I- s^Aa — spc^^X^ — 2 s^c^ReA5 -1- ReAg -|- c^(c^ — 2s^) ReA 7 , 

gHH+H _ 2s^c^ImA5 — s^ImAe — c^ImAr, (A.l) 


^111 = cgXi + -s^aReAe, 

3 

9 \i 2 — -^34 + sg ReAs + -cp ReAg , 

3 

9i22 — C/3 A 34 -|- C /3 ReAs + -;^sp ReAr , 

9222 = sp^2 + 2 ReAr , 

fi '113 = -spCplmX^ - ^ (1-F 2cJ) ImAe , 


9 I 2 Z 

9^23 


9i33 


9 I 33 


9 I 33 


-2ImA5 - spcplui (Ae -h A7), 

-spcp IniAs - ^ (1 -F 2sp) IniA7 , 

s'pCpXi + A34 — C/3(l-|-s^) ReAs + ~2^p^‘p^^^7 , 

spc‘pX 2 + A34 — 5 / 3 ( 1 -f-c^) ReAs + T 2 ~ ’ 

spCplmX^ - ^sJlmAe - ^ cj ImA7 , (A.2) 


with A 34 = I (A 3 -|- A 4 ) and sp (cp) = sin/3 (cos/3). Onr analytic expressions in (|A.1|) 
and (|A. 2 |) agree well with those presented in [p7[|.P| 


In the remainder of the Appendix, we present new analytic resnlts for the proper 
qnadrilinear self-conplings and ■ These are given by 


9uii = 


9ui2 = 

2 ReAe , 

9u22 = 

iA34 + 

— ReAs 
2 

9t222 = 

2 

9t222 






^Recently, another study of Higgs self-couplings appeared |Q. The authors expressed their analytic 
results in terms of CP-conserving mixing angles and Higgs-boson masses, thus rendering a direct comparison 
with our expressions very difficult. 
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£/m3 


9t223 


9u33 


9^233 


92233 


9u33 


92^33 


93^33 


--c^ImAe, gn 23 = -cglmAs - - ImAg , 

-S/jImAs - ^C/jImAr, 5(3^3 = -^s^lmAr, 

- Ai + - A34 ~ 2 ReAs — - s^C/3 ReAg, 

—2s/ 3C/3 ReAs 2 ^9 ReAg + - ReA7, 

111 1 

- A2 + - A34 ~ 2 ^‘9 ReAs “ 2 ^ 9^9 ReAy , 

S/3cJlmAs - ImAe - ^cJlmAy, 

sJc/jImAs - ^sJlmAe - ImAy, 

1 ^ 4 ^+ 

4 


(A.3) 


where has been presented in ([^.341), and 


^2HH+H~ 

9n 

— Ai + - A 3 — spCp 

^2HH+H- 

9i2 

= —'S/ 3 C/ 3 A 4 — 2s/3C/3ReAs 

^2HH+H- 

922 

= A 2 + - A 3 — S/ 3 C /3 

^2HH+H- 

9i3 

= 2sp(?p ImAs — s'^pCp ImAf 

^2HH+H- 

923 

= 2s'pCp ImAs ~ Sp ImAg • 

2HH+H- 

933 



(A.4) 


Notice that the proper self-couplings and gf^ may be also expressed in terms of 

gfHH+H- follows: 


9HH^h- 


„ 2HH+H- , 

C/3 9i3 + 


9iii — 9iiii + 'S/3 5'iii2) 

^112 = 3c^^m2 + 

5'122 ~ ^^g9l222 R gii22 I 

„3H _ A AH I AH 

9222 ~ ^*/9 92222 ' C /3 (/i 222 ) 

9u3 = 3c/3C/m3 + S/3^m3> 

9i23 ~ ^C/3 (71123 + ^'^/3 5'1223 ’ 

_ o _|_ 

i/223 “ i/2223 i/12235 


(A.5) 
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9 i 33 ~ (71133 + ^0 9 i 233 > 

fi'233 ~ (72233 + fi'1233 5 

fi'333 ~ ^0 9 i 333 + ^0 92333 ' (^-6) 


The above relations are useful for checking the self-consistency of all the analytic expressions 
pertinent to the trilinear and quadrilinear Higgs self-couplings. 
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2 Experimental Simulations used in our Analysis 


We summarize here the results from previous simulations of Higgs production at the LHC, 
upon which we have based our analysis of MSSM neutral-Higgs production in the presence 
of CP violation. In each Table, we provide the minimum ratio R, relative to the Standard 
Model product of cross section times branching ratio, for which a signal of the stated signif¬ 
icance could be obtained with the quoted luminosity. We have provided more information 
regarding the expected number of signal S and background B events for those cases when 
S/B~ 1. 


Mh (GeV) 

90 

100 

110 

120 

130 

R3.(5 fb-i) 

0.63 

0.71 

0.88 

1.07 

1.43 


Table 1: R values for 3-a evidence for a Higgs boson in the channel W/ZH{^ bb) at the 
Tevatron with 5 fb~^ of accumulated data ]^/ . These results are based on combining the 
CDF and D0 data. 


Mh (GeV) 

S (100 fb-^) 

B (100 fb-^) 

^Gauss 

Rs. (100 fb-^) 

100 

147 

223 

9.84 

0.51 

105 

123 

191 

8.92 

0.56 

110 

100 

160 

7.91 

0.63 

115 

98 

163 

7.68 

0.65 

120 

96 

167 

7.43 

0.67 

125 

86 

157 

6.89 

0.73 

130 

77 

148 

6.33 

0.79 


Table 2: The expected signal, background, Gaussian significance and the R value for 5-a 
discovery of a Higgs boson in the channel tiH{-^ bb) at the LHC with 100 fb~^ of accumu¬ 
lated data These results are for the CMS experiment only. The Gaussian significance 
is based on a K factor of unity for the signal j^J . 
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Mh (GeV) 

S (30 fb-i) 

B (30 fb-i) 

^Gauss 

Rscr 

110 

11.1 

3.9 

4.1 

1.22 

120 

10.4 

1.4 

5.2 

0.96 

130 

8.6 

0.9 

5.0 

1.00 

140 

5.8 

0.7 

3.9 

1.28 

150 

3.0 

0.6 

2.3 

2.17 


Table 3: The expected signal, background, Gaussian significance and the R value for 5-a 
discovery of a Higgs boson in the channel WW —>• H{-^ t^t~) at the LHC with 30 fb~^ of 
accumulated data j^/ . A detector-level analysis based on a fast simulation was presented 
in |g^/ . The results there for the combination of one hadronic and one leptonic decay 
of the T^T~ pair are slightly less significant than the original analysis, but they have not 
been combined with the all leptonic decays. Therefore, we have chosen to use the original 
numbers, even though they are based on a parton-level analysis. 


Mh (GeV) 

80 

90 

100 

110 

120 

130 

150 

Rs. (100 fb-i) 

1.02 

0.74 

0.60 

0.49 

0.42 

0.4 

0.66 


Table 4: The expected signal, background, Gaussian significance and the R value for 5-a 
discovery of a Higgs boson in the channel gg —> H{^ 77) at the LHC with 100 fb~^ of 
accumulated data These apply to CMS only. 
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